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ABSTRACT 

A convincing explanation for the observed pulsar large peculiar velocities is still missing. We 
argue that any viable particle physics solution would most likely involve the resonant production 
of a non-interacting neutrino v s of mass m Va ~ 20-50 keV. We propose a model where anisotropic 
magnetic field configurations strongly bias the resonant spin flavour precession of tau antincu- 
trinos into v s . For internal magnetic fields B lnt > 10 15 G a v T -v s transition magnetic moment 
of the order of 10 -12 Bohr magnetons is required. The asymmetric emission of v s from the core 
can produce sizeable natal kicks and account for recoil velocities of several hundred kilometers 
per second. 

Subject headings: elementary particles — pulsars: general — stars: neutron — supernovae: general 



1. INTRODUCTION 

Measurements of the pulsars proper motion 
have shown that the pulsar population is char- 
acterized by peculiar velocities much higher than 
those of all other stellar populations (Gunn & 
Ostriker 1970; Lyne & Lorimer 1994; Hansen & 
Phinney 1997; Cordes & Chernoff 1997; Lorimer, 
Bailes, & Harrison 1997). While the precise form 
of the velocity distribution is uncertain, recent 
studies indicate that it must extend to ~ 1000 km 
s _1 , with a mean of ~ 500 km s _1 (Lorimer et 
al. 1997). Since the velocities of pulsars' likely 
progenitors, the massive stars, are among the low- 
est in the galaxy (~ 10-40 km s" 1 ) it is natu- 
ral to think that the mechanism responsible for 
their acceleration can be related either to super- 
nova (SN) core collapse and explosion, or to the 
post-explosion cooling phase when, within a few 

1 Permanent address: Departamento de Fi'sica, Univer- 
sidad de Antioquia, A. A. 1226, Medellm, Colombia 



seconds, a huge amount of energy is emitted from 
the proto-neutron star (PNS) through the radia- 
tion of ~ 10 58 neutrinos of all flavours. 

To date, no convincing physical process that 
could account for pulsar acceleration has been 
identified (a review of different kick mechanisms 
has been recently given by Lai (1999)). The as- 
trophysical solutions that have been proposed in- 
clude, in broad terms, the possibility of symmetric 
and of asymmetric SN explosions. In the first case, 
large recoil velocities are achieved in response to 
the mass loss in the (symmetric) explosion of one 
component of a close binary system (Gott, Gunn, 
& Ostriker 1970; Iben & Tutukov 1996, 1997). It 
was shown that the remnant can achieve velocities 
comparable to the orbital velocity of the SN pro- 
genitor. However, detailed Monte Carlo simula- 
tions of binaries' evolution through the SN phase 
fail to reproduce the observed velocity distribu- 
tions of pulsars and hint at asymmetric explo- 
sions (Lorimer et al. 1997; Dewey & Cordes 1987; 



van den Heuvel & van Paradijs 1997; Hughes & 
Bailes 1999). In the second case, asymmetrical 
stellar collapse and explosion occurs because of 
pre-explosion instabilities and core density distor- 
tions, imparting to the residue an intrinsic kick 
(Janka & Mueller 1994; Burrows & Hayes 1995, 
1996). The momentum scale of the mass ejecta in 
the SN explosion is of the order of w 6x 10 53 erg/c, 
corresponding to about ten solar masses ejected 
with velocities 10 4 km s _1 . Then, to acceler- 
ate the remnant of about 1.5 M Q to the largest 
observed velocities, it is necessary to generate a 
momentum asymmetry in the ejecta of the order 
of a few per cent. The asymmetric explosion sce- 
nario remains unsatisfactory in the sense that ini- 
tial anisotropies in the collapsing core sufficiently 
large to generate asymmetries of this size have to 
be input artificially (Burrows & Hayes 1996) and, 
from the point of view of present (2-D) hydrody- 
namic calculations, are difficult to justify (Janka 
& Mueller 1994). 

More recently it has been suggested that a com- 
pletely different class of mechanisms could be re- 
sponsible for imparting a natal kick to the neu- 
tron star. Virtually all the gravitational binding 
energy released during the collapse of the progen- 
itor iron core (w 3 x 10 53 ergs) is emitted during 
the first 10 seconds after core bounce in the form 
of neutrinos and antineutrinos. The momentum 
in neutrinos is therefore comparable to the mo- 
mentum in the ejected matter. Hence, also in this 
case a few per cent asymmetry in the neutrino 
emission could suffice to account for the largest 
velocities observed. The particle physics explana- 
tions for the pulsar acceleration mechanisms can 
be classified according to the underlying mecha- 
nism at the origin of the momentum asymmetry 
in the neutrino emission. The asymmetry could 
be in the number of neutrinos, or in the neutrino 
energy (since the number of neutrinos is not con- 
served when neutrinos stream out of the core, this 
only refers to the microscopic cause of the asym- 
metry, and does not describe the macroscopic ef- 
fect). In both cases, the presence of a strong mag- 
netic field B — 10 15 -10 16 G provides the asym- 
metric initial condition for the system. In the first 
case, parity violation in clectroweak interactions 
was thought to be at the origin of asymmetric 
neutrino emission from strongly magnetized PNS. 
The cumulative effect of multiple neutrino scat- 



terings off the slightly polarized nucleons would 
result in a sizeable anisotropy in the neutrino mo- 
mentum (Horowitz & Piekarewicz 1998; Horowitz 
& Li 1998; Lai & Quian 1998a). More detailed 
studies showed that this effect was largely overesti- 
mated (Kusenko, Segre, & Vilenkin 1998; Arras & 
Lai 1999). When neutrinos are in local thermody- 
namic equilibrium (as they are to a very good ap- 
proximation in the interior of a PNS) detailed bal- 
ance ensures that no cumulative effect from mul- 
tiple scatterings can build up. Local effects in the 
PNS atmosphere, as for example the deviations 
from thermal equilibrium that occurr close to the 
neutrinosphere, can generate asymmetries in the 
neutrino emission. However, to date detailed anal- 
yses of these effects have failed to produce kicks as 
large as 1000 km s" 1 (Arras & Lai 1999). The sec- 
ond case relies on the assumption that resonant 
transitions (oscillations (Kusenko & Segre 1996, 
1997; Grasso, Nunokawa, & Valle 1998) or spin- 
flavour precession (Akhmedov, Lanza, & Sciama 
1997)) between different neutrino flavours occur 
inside the PNS. The resonance surface for the 
conversion becomes the effective surface of last 
scattering (the neutrinosphere) for the neutrino 
flavour with the largest mean free path. Since a 
strong magnetic field can slightly distort the reso- 
nance surface, neutrinos are emitted from regions 
with different temperatures, which yields an asym- 
metry in the overall momentum. Also in this case, 
the real effect was largely overestimated. This 
is because the variation of the temperature over 
the effective neutrinosphere cannot be directly re- 
lated to an anisotropy in the neutrino energy. The 
energy flow from the star is governed by the in- 
ner core emission of neutrinos; in first approxi- 
mation, local processes in the PNS atmosphere, 
where the resonant conversion takes place, cannot 
modify the (isotropic) neutrino emission from the 
core (Janka & Raffelt 1999). 

In conclusion, both the macroscopic parity vi- 
olation and the resonant neutrino conversion sce- 
nario, which for a while seemed to provide elegant 
particle physics mechanisms to account for pul- 
sars' acceleration, did not survive more detailed 
analysis. However, the study of these attempts 
and of the reasons why they fail to produce natal 
kicks of the right magnitude can guide us to infer 
some of the conditions that must be satisfied by a 
viable particle physics solution. 
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• The neutrino observations from SN1987A, 
and in particular the duration of the v e 
burst, support the theoretical expectation 
that active neutrinos are copiously pro- 
duced inside the PNS core and that they 
are the main agents of the energy emis- 
sion. Theoretical studies indicate that, to a 
good approximation, the total energy carried 
away by neutrinos is equipartitioned among 
the six neutrino and antineutrino flavours 
(Janka 1995). This implies that conversions 
between the active flavours would simply 
result in exchanging the respective spectra, 
and cannot affect the main characteristics of 
energy emission (Janka & Raffelt 1999). We 
learn from this that if neutrino conversion is 
responsible for the kicks, then the conversion 
must be into some new particle that does not 
interact with matter in the same way than 
the standard neutrinos: most probably a 
very weakly interacting or sterile neutrino 
v s . 

• In thermal equilibrium, the microscopic pe- 
culiarities of neutrino scattering processes 
cannot build up any macroscopic asymme- 
try. This applies also to possible neutrino 
scattering into new exotic states. This hints 
at some macroscopic condition anisotropi- 
cally realized inside the PNS. If this has to 
influence the neutrino conversion rates, we 
are naturally led to thinking at some asym- 
metrically located resonance, or at local vi- 
olation of the adiabaticity conditions for the 
conversion. 

• An anisotropic resonant conversion of an ac- 
tive neutrino into an exotic state is still not 
sufficient to generate a momentum asymme- 
try. Close to the resonance region the parent 
v and the v s would both have an asymmetric 
momentum distribution, but their asymme- 
tries would compensate in the total momen- 
tum. However, a non-interacting v s would 
freely escape from the PNS, thus preserving 
its original asymmetry, while the active neu- 
trino keeps interacting with the background 
particles. If v has enough time to recover 
a symmetric momentum distribution before 
escaping from the star, then the v s asymme- 
try will be left unbalanced. This suggests 



that the conversion should occur far inside 
the v neutrinosphere, and most likely in the 
high-density regions close to the core. 

The picture emerging from these considerations 
is a resonant conversion of an active flavour into 
a non-interacting neutrino, at the large densities 
and temperatures typical of the PNS regions close 
to the core. In turn, matching the resonance con- 
dition hints at a sterile neutrino mass « 20-50 
keV. In order to implement the previous scheme 
it is important to identify all the possible macro- 
scopic asymmetries of the system. As we will dis- 
cuss in the next section, there are observational 
indications that, in some cases, pulsar magnetic 
fields are better described by an off-centred dipole, 
rather than simply by an oblique dipole with the 
centre of the magnetic axis coinciding with the 
centre of the star. Clearly this implies anisotropic 
conditions inside the PNS core, since the mag- 
netic field strength can easily vary by more than 
one order of magnitude within regions of equal 
matter density. Then it is not difficult to envis- 
age how in such an anisotropic background the 
Mikheyev-Smirnov-Wolfenstein (MSW) mecha- 
nism for matter-enhanced neutrino oscillations 
(Wolfestein 1978; Mikheev & Smirnov 1985) or 
resonant spin-flavour precession (RSFP) (Akhme- 
dov 1988a, 1988b; Lim & Marciano 1988) could 
result in large asymmetries in the rate of neutrino 
conversion between opposite hemispheres. More- 
over, if the conversion involves non-interacting 
neutrinos, they will freely escape from the core, 
yielding a large momentum asymmetry. We stress 
that the solution we are proposing here is intrin- 
sically different from previously studied mecha- 
nisms that exploit resonant neutrino conversion 
and emission from a slightly deformed resonance 
surface (Kusenko & Segre 1996, 1997; Grasso et 
al. 1998; Akhmcdov ct al. 1997). In our case, it is 
the conversion probability that is not isotropically 
distributed inside the (spherical) region where the 
resonant condition is satisfied. 

In Section 2 we will review some of the argu- 
ments in support of our basic assumption that the 
magnetic field in the PNS interior can be charac- 
terized by sizeable anisotropics. In Section 3 we 
will briefly review the physics of RSFP. For defi- 
nitencss, we will assume that tau antineutrinos can 
be resonantly converted into sterile neutrinos v s . 
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We will show that for acceptable values of the PNS 
magnetic field and of the neutrino transition mag- 
netic moment \x v the adiabaticity conditions for an 
efficient conversion are satisfied in regions asym- 
metrically displaced from the star centre. A MSW 
solution to the problem can also be implemented 
straightforwardly since the asymmetrically polar- 
ized background results in an anisotropic neutrino 
potential (Nunokawa et al. 1997). However, this 
requires more fine tuning in the choice of the rele- 
vant parameters and we will only comment briefly 
on this possibility. The details of the PNS model 
we have used and of neutrino diffusion from the 
core will be described in Section 4. Finally in Sec- 
tion 5 we will present our results and draw the 
conclusions. 

2. EVIDENCE FOR ANISOTROPIC MAG- 
NETIC FIELDS 

The evidence in support of asymmetric mag- 
netic field topologies for variable magnetic stars 
and pulsars was discussed long ago by Harrison 
and Tademaru (Harrison & Tademaru 1975; Tade- 
maru & Harrison 1975; Tademaru 1976). They ar- 
gued that pulsars' magnetic fields can be well mod- 
eled by an off-centred magnetic dipole, displaced 
from the centre of the star and oriented obliquely 
with respect to the rotation axis in an arbitrary 
direction. They also showed that such a field con- 
figuration radiates asymmetrically low-frequency 
electromagnetic radiation, and they tried to ex- 
plain through this mechanism 'post-natal' pulsar 
acceleration at the expense of the star rotational 
energy. However, velocities much in excess of 100 
km s _1 cannot be explained in this way since for 
typical rotational periods of the order of 10 ms 
the rotational energy does not exceed by much 
the total kinetic energy of the star. Moreover, the 
lack of any observed correlation between the mag- 
netic field and the transverse velocity, or between 
the direction of the motion and the magnetic axis 
(Lorimer et al. 1997; Deshpande, Ramachandra, & 
Radhakrishnan 1999) has virtually ruled out the 
post-natal acceleration model. Asymmetric mag- 
netic field configurations can also affect v e and v e 
absorption cross sections on neutrons and protons. 
Lai and Qian (1998b) studied the possible asym- 
metry in the neutrino emission that can be gen- 
erated in this way and showed that the effect is 
again too small to produce appreciable kicks. 



However, an anisotropic magnetic field can still 
be responsible for pulsar acceleration by producing 
sizeable asymmetries in the neutrino emission dur- 
ing the early PNS cooling phase. In principle, the 
mechanism we want to propose can be effective in- 
dependent of the particular type of magnetic field 
model, as long as it predicts sizeable anisotropics 
in the field strength. However, for definiteness we 
will concentrate on the simple off-centred and de- 
centred magnetic dipole models (Harrison & Tade- 
maru 1975; Tademaru & Harrison 1975; Tademaru 
1976) also because they are supported by a few 
observational evidences. Decentred dipole models 
were first introduced in the attempt to improve 
the representation of the fields of periodic mag- 
netic variable stars (Landstreet 1970). It was later 
suggested that these models could also account for 
special features in the electromagnetic emission of 
the so-called intcrpulse pulsars, which are char- 
acterized by secondary interpulses well separated 
from the main pulse, and of an intensity up to 
two orders of magnitude smaller (Hankins 1986). 
A displacement of the magnetic dipole along the 
magnetic axis of an amount 5 in units of stellar 
radius would produce unequal surface field inten- 
sities in the ratio (1 + S) 3 /(l — 5) 3 , providing a 
simple explanation for the relative weakness of 
the secondary interpulse. The decentre parame- 
ter (5 inferred from observational studies of vari- 
able magnetic stars and interpulse pulsars ranges 
between 0.1 and 0.6 (Harrison & Tademaru 1975; 
Tademaru & Harrison 1975; Tademaru 1976). In 
some cases surface fields of opposite polarities, 
separated by much less than 180° have been ob- 
served (Hankins 1986). Asymmetries of this kind 
are not explained by the simple decentred model, 
and suggest that the dipole is actually off- centred. 
Namely, the magnetic dipole \x = (/i z , /i p , /i^) (in 
cylindrical coordinates) is arbitrarily oriented and 
displaced from the centre of the star by an amount 
5 Z along the spin axis and 6 P in the radial direc- 
tion. Since the dipole axis does not intersect the 
spin axis, the angular separation between the main 
and the secondary pulses is easily accounted for. 

The study of interpulse pulsars is statistically 
limited by the small sample of a few per cent of 
all pulsar population. This is likely to be due to 
the fact that both the angle of inclination between 
the rotation axis and the line of sight, and the an- 
gle between the spin and the magnetic axis, must 
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be close to 90° to allow for their detection. Still, 
the few direct observations available today suggest 
that it is not unreasonable to assume that most 
pulsars have oblique and off-centred dipolc fields, 
or possibly some more complicated field configura- 
tion implying sizeable anisotropics in the magnetic 
field topology. 

Since off-centred dipoles radiate energy at a dif- 
ferent rate with respect to centred dipoles, one 
might wonder if off-centred models could be tested 
by measuring the pulsar rate of energy loss. In 
simple cases, the loss of angular velocity w of a 
pulsar can be described by the law 

Cb = -k(O n , (1) 

where n , called the braking index, can be directly 
inferred from the instantaneous rate of change of 
frequency through the relation n = Colo j Co 2 . A 
magnetic dipole /x located on the spin axis radi- 
ates electromagnetic energy at the expense of ro- 
tational energy at a rate 

dE^ 1 djlco 2 ) _ 2y 2 4 
dt 2 dt 3c 3 ^ ' U 

where I is the neutron star moment of inertia. For 
constant / this yields n — 3. For gravitational 
quadrupole radiation, denoting by M p the com- 
ponent of the mass-quadrupole orthogonal to the 
spin axis, dE M /dt — (GnM 2 /45c 5 ) lo 6 is found 
(Ostriker & Gunn 1969), and the braking index is 
n = 5. For a magnetic dipole displaced from the 
spin axis by an amount 6 p R , the spin-down law 
reads 

Co = -ku 3 (l + aw 2 ) (3) 

where k = 2/^/3/c 3 and a = 2y 2 z {5 p R) 2 /by^c 2 . 
This yields a braking index n — Colo /Co 2 ~ 3 + 
2aio 2 > 3. Unfortunately, for typical values of the 
relevant parameters (P <~ 10 ms, R ~ 10 km, etc.) 
one obtains au; 2 < 10~ 5 , and thus the deviation 
from the centred dipole spin-down law is by far 
too small to be measurable. 

3. RESONANT SPIN-FLAVOUR PRE- 
CESSION 

The mechanism we want to explore to account 
for pulsar acceleration is the RSFP of an active 
neutrino into a sterile neutrino, driven by the 
anisotropic magnetic field of a decentred dipole 



configuration. For simplicity, we will not analyze 
here the more realistic case of an off-centred mag- 
netic dipole or of more complicated configurations, 
since these generalizations do not change our main 
results. We just note that, in the off-centred case, 
sterile neutrino emission could contribute not only 
to the acceleration, but also to a spin-up of the 
PNS. As has recently been suggested (Spruit & 
Phinney 1998) the rotating core of the pulsar pro- 
genitor might have just a tiny fraction of the an- 
gular momentum required to explain the observed 
rotation velocities of pulsars. For this reason it 
was conjectured that the same physical process 
that kicks the PNS at birth could be responsible 
also for pulsars' fast rotations (Spruit & Phinney 
1998; Cowsik 1998). Thus, it might be interesting 
to study if a more general magnetic field configu- 
ration, in conjunction with neutrino RSFP, could 
also account for the observed pulsars periods. 

We start by assuming that v T — ► v s RSFP can 
occur at neutron star core densities and magnetic 
field strengths. Since v s can freely escape from the 
core, we need to ensure that the neutrino burst 
duration ~ 10 s will not be shortened too much. 
Assuming equal luminosity for all neutrino and an- 
tineutrino flavours (Janka 1995), even a 100% effi- 
cient conversion of v T into v s would affect only 
1/6 ~ 17% of the energy in the game. Since 
there are resonance regions where the adiabatic- 
ity conditions are not satisfied and no conversion 
occurs, the overall efficiency for v T — > v s conver- 
sion is in fact lower, meaning that the amount 
of energy carried away by v s is small enough for 
the burst duration not to be drastically changed. 
On the other hand a 10%-20% asymmetry in the 
active-sterile conversion efficiency between oppo- 
site hemispheres would result in a few per cent 
asymmetry in the total momentum, as is needed 
to explain the largest velocities observed. An im- 
portant point in this discussion is to ensure that 
transitions that involve other active flavours, such 
as v e . p — ► v T — > v s , will not result in an en- 
ergy 'siphon' effect that would cool the PNS too 
rapidly. Indeed, in the light of the present ex- 
perimental hints for non-vanishing neutrino mix- 
ings, accounting for these reactions is mandatory. 
The rates for neutrino flavour conversion in a neu- 
tron star core were studied by Raffelt and Sigl 
(1997) and by Hannestad et al. (1999). For v e 
the large charged current (CC) refractive effects 
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due to the background electrons strongly suppress 
any in-matter mixing angle, so that conversion to 
other neutrino species cannot occur on the time 
scale of neutrino diffusion (Raffelt & Sigl 1993). 
Also for Vy, <-» v T the conversion time scale safely 
exceeds the neutrino diffusion time (Hannestad 
et al. 1999). This is due to the presence of a 
small amount of muons in the hot superdense core, 
which implies additional CC contributions to the 
Vy index of refraction, as well as to the second- 
order difference in the Vy-v T neutral current (NC) 
potentials SV ~ G^m 2 A^v (where Nn is the den- 
sity of nucleons) (Botella, Lim, & Marciano 1987). 
In conclusion, the large differences in the neutri- 
nos indices of refraction inside the PNS guarantee 
that e, \x and r lepton numbers are separately con- 
served, and thus the conversion into v s will affect 
just one neutrino flavour. 

The condition for a RSFP of v T into sterile neu- 
trinos v s reads 

T ,,_ x G F N n ^ f Am 2 

f=n,p,e 

The effective potential V(D T ) felt by a v T prop- 
agating in the hot, superdense and slightly po- 
larized PNS matter has two contributions: the 
first one is due to the coherent vector interac- 
tion with the background and is proportional to 
the number density of neutrons N n . The second 
one, which is due to coherent axial-vector inter- 
actions, is proportional to the average polariza- 
tion (Aj|) of the / = n,p,e background fermion 
parallel to the direction of the neutrino motion 
(Nunokawa et al. 1997; Bergmann, Grossman, & 
Nardi 1999). The proportionality factor kf ~ 
id A Gp/VtyNf depends on the/-fermion number 
density Nf and on its axial- vector coupling g A . 
In the right-hand side of equation (4), Am 2 de- 
notes the square mass difference between the tau 
and the sterile neutrino masses. In the PNS core 
V{v T ) w 4 x lCT 6 /Ji4MeV, where p u = p/10 14 g 
cm~ 3 . For neutrino thermal energies of the order 
E v ~ 100-200 Me V the resonance condition can 
be satisfied if Am 2 - (20-50 keV) 2 . This value 
is much larger than the cosmological limit to„ t 
(< 10 2 cV), and in the following we will therefore 
assume Am 2 = m 2 — m 2 ~ m 2 . 

is s i/ T 1/ s 

The average background polarization con- 
tributing to the second term in V(P T ) grows lin- 



early with the magnetic field, which in a neu- 
tron star can reach extremely high values, up 
to - 10 15 G (Thomson & Duncan 1992, 1993, 
1995, 1996). However, during the early cooling 
phase when most of the neutrinos are emitted, 
the temperature is also large, the electrons are 
relativistic and degenerate, and as a result the in- 
duced polarizations are strongly suppressed, down 
to (Vp> < 1CT 3 x B 15 and (A e > < 10~ 2 x B 15 
(Bergmann et al. 1999), where B15 is the mag- 
netic field strength in units of 10 15 G. Thus in 
general K/(Aj|) <C GfN u /\/2 and, in most cases, 
neglecting the polarization term is well justified. 

The resonance condition for the (hclicity- 
conserving) MSW effect (p T <-> v s ) has the same 
form that equation (4) except for the fact that the 
right-hand side is multiplied by cos 29y , where 9y 
is the vacuum mixing angle. Because the magnetic 
dipolc configuration is asymmetric, in different re- 
gions of a core shell of given neutron density N n 
the magnetic field and the average polarization 
(A-^) can be vastly different. Then it could be pos- 
sible that the resonance condition is satisfied with 
sufficient accuracy in just one hemisphere. This 
can be more easily achieved for resonant D e <-> v s 
oscillations, since in this case the neutron number 
density in the first term in equation (4) is replaced 
by iV n — 2N e , and in some region of the core and 
at some stage of the PNS evolution this combina- 
tion can approach values close to zero. However, 
local violation of the adiabaticity condition for a 
RSFP conversion (see below) appears to be more 
natural than an asymmetric MSW resonant con- 
version, in the sense that it can be satisfied for a 
larger range of the relevant parameters, and is also 
more stable with respect to core evolution during 
the deleptonization and cooling phases. Therefore 
in the following we will concentrate on the RSFP 
conversion mechanism. 

Once the condition in equation (4) is satisfied, 
the probability P^ T ^ lJs that a spin-flavour neu- 
trino transition will occur depends on its degree 
of adiabaticity. To a good accuracy 

Pv T ^u 3 ^ 1 - exp 7) , (5) 

which is sizeable when 7 > 1. Denoting by 
£ p = \(l/p)(dp/dr)\^ 1 the characteristic length 
over which the density varies significantly, and 
with £p r cs its value at the resonance, the adia- 
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baticity parameter 7 can be written as (Akhmedov 
1997; Akhmedov et al. 1997) 



7 ~ 1.04 



1 \ /2.6 x 10 14 gcnT 



l-Y P 



fjb v B 



Xros 



10- 12 ms 3 x 10 15 G/ 10km 



(6) 



where Y e w 0.3 is the number of electrons per 
baryon in the PNS core (Burrows & Lattimer 
1986), fiB is the Bohr magneton and B r f a is the 
value at resonance of the magnetic field component 
orthogonal to the direction of neutrino propaga- 
tion. Let us now assume that the magnetic dipolc 
is displaced from the star centre by an amount S 
in units of the PNS radius R , and let us define 
5 = 5- (R/R Tes ) , where R ICS is the location of the 
resonance layer. For reasonable values 5 w 0.2- 
0.5 the ratio of field intensities between the north 
(N) and south (S) magnetic poles of the resonance 
shell B T s es /B™ s ~ (l-S) 3 /(l+5) 3 falls in the range 
10 -1 -10 -2 . Then from equation (6) we see that 
the condition 7jv > 1 (Pv T ^ Va ~ 1) and 75 < 1 
(Pp T< _>i, s ~ 0) that triggers the anisotropic neu- 
trino emission can be realized in a natural way. 
In Fig. 1 we depict the variation of the adia- 
baticity parameter 7 inside the resonance layer 
as a function of the angular distance 6 from the 
magnetic dipole axis. The resonance is located 
at i? rC s = 1-5 R c , where R c denotes the core ra- 
dius, and the maximum value of the magnetic field 
strength at resonance is B = 4 x 10 15 G. For small 
values of the decentred parameter S ~ 0.1 the 
adiabaticity condition is matched in both hemi- 
spheres and no large asymmetry can be expected. 
However, for S > 0.2 we get 7 < 1 in the whole 
hemisphere 8 > ir/2 and a sizeable asymmetry 
can be produced. It is also apparent that if S be- 
comes too large, the region where the adiabaticity 
condition is satisfied shrinks down to a small cone 
(e.g. for S > 0.5 there is good adiabaticity only 
when 6 < 7r/6). This heavily reduces the conver- 
sion efficiency and results in a suppression of the 
asymmetry. It also suggests that no simple cor- 
relation between the value of S and the overall v s 
momentum asymmetry can be expected. 

Let us now comment briefly on the possible 
values of the relevant physical quantities appear- 
ing in equation (6). Magnetic field strengths 
> 10 15 G arc not unreasonable in the interior of 



a new-born neutron star. Indeed, there is ob- 
servational evidence for highly magnetized young 
pulsars (magnetars) with dipole surface fields as 
large as 8 x 10 14 G (Hurley 1999). On the other 
hand, theoretical studies indicate that internal 
dipole fields as large as B int — (5-10) xl0 15 G 
can be formed during the first few seconds after 
gravitational collapse, when the strong convective 
motions coupled with rapid rotation can produce 
an efficient dynamo action (Thomson & Duncan 
1992, 1993, 1995, 1996). As regards the value of 
[i v , it is well known that the RSFP conversion 
mechanism requires rather large transition mag- 
netic moments. For densities close to nuclear den- 
sity (w 2.6 x P14) and acceptable values of the 
magnetic field, the adiabaticity condition in equa- 
tion (6) points towards magnetic moments of the 
order of 10~ 12 ^b. For v T this is several orders 
of magnitude below the laboratory limits (Groom 
et al. 1998). In the early Universe, v s will attain 
thermal equilibrium only for [i v > 6 x 10~ n ^-B 
(Elmfors et al. 1997) and hence the nucleosynthe- 
sis constraints can be easily satisfied. The astro- 
physical limit [i v < 3 x 10~ 12 implied by stel- 
lar energy-loss arguments (Raffelt 1990) does not 
apply in this case, since the final-state sterile neu- 
trino with m s » 5 keV is too heavy to be produced 
inside red-giants or white-dwarfs. Assuming, as is 
reasonable, that the main v s decay mode is radia- 
tive v s — > P T 7 , the limits on 7-ray fluence imme- 
diately after the arrival of the neutrino burst from 
the SN1987A can be used to set strong constraints 
on fi v . Under the assumption that the decaying 
neutrinos are carrying away about 1/3 of the total 
energy, the limit \x v < 1.6 x 10~ 14 /Z-B (10keV/m„ s ) 
was derived (Oberauer et al. 1993; Raffelt 1996, 
p. 474). However, this limit cannot be straightfor- 
wardly applied to our case. As we will see below, 
depending on the overall conversion efficiency, the 
v s can easily carry an overall energy more than 
one order of magnitude smaller than what was 
assumed in order to derive this limit. It is even 
conceivable that no conversion at all could have 
occurred for SN1987A, if the magnetic field was 
too weak to satisfy the adiabaticity condition (of 
course this would also imply no intrinsic kick for 
the SN1987A residue, which unfortunately has not 
been detected). In addition to this, in our frame- 
work the v s emission is strongly anisotropic, and 
beaming effects could have drastically reduced the 
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flux along the line of sigh of the Earth. We con- 
clude that the SN1987A limit cannot exclude val- 
ues of p, v of the order 10~ 12 /i_b- However, it is 
interesting to note that a clear signature of the 
mechanism we are proposing would be a large flux 
of 7 rays with energies of several tens of MeV from 
the next galactic SN. Since the v a are emitted from 
regions close to the hot inner core, this would also 
be accompanied by an anomalous spectral com- 
ponent of v T from v s decays with a temperature 
about one order of magnitude larger than the v r 
and Vy, spectra. 

If f ew x10 _12 ^b the sterile neutrinos 

would be directly produced through helicity flip- 
ping scattering off electrons and protons (Barbi- 
eri & Mohapatra 1989; Ayala, D'Olivo, & Torres 
2000). The main concern here is not so much 
the change in the rate of energy loss, since this 
would affect only one flavour, but rather the fact 
that spin-flips due to scattering processes would 
then occur within the entire core volume with an 
(almost) isotropic distribution. If the time scale 
for these processes is much shorter than the v T 
diffusion time to the resonant region, almost all 
the v T will convert into v s inside the core. How- 
ever, in crossing the resonance layer, the v s will be 
anisotropically reconverted into interacting states. 
This would result in approximately the same mo- 
mentum asymmetry, but in the opposite direction. 
If the time scale for the helicity flipping scatterings 
is much larger than the v T diffusion time to the res- 
onance (but of course shorter than their diffusion 
time to the neutrinosphcre) the v s asymmetry gen- 
erated at the resonance surface will be preserved. 
This is because the unconverted D T streaming out 
from the resonance layer will quickly recover an 
isotropic momentum distribution, so that the later 
emission of v s from helicity flipping scatterings 
will be essentially symmetrical. Still, allowing for 
these possibilities would unnecessarily complicate 
our analysis. Therefore, to ensure that helicity 
flipping scatterings are safely suppressed and will 
not interfere with resonant conversions, we will as- 
sume that the limit \i v < few xlO~ 12 yU_B (Barbieri 
& Mohapatra 1989; Ayala et al. 2000) is satisfied. 

4. THE PROTO-NEUTRON STAR MODEL 

In this section we describe the model that was 
used to simulate the PNS physical conditions dur- 



ing the first few seconds after core bounce. Our 
PNS model is three dimensional but static, so that 
the results we will obtain should be understood 
MS M "proof-of-principle" calculation of the possi- 
ble size of the momentum asymmetry. We iden- 
tify four different regions that are of major im- 
portance for the conversion process and for the 
neutrino emission: 

1) Core- atmosphere interface (r = R c ). Neu- 
trinos are mainly produced within the hot 
and dense core, where the density varies 
slowly with the radius. From the central re- 
gions with supernuclear density p a ~ 8 x pi4 
the density decreases down to nuclear den- 
sity p c ~ 2.6 x P14. We define the core- 
atmosphere interface radius R c through 
p{R c ) — p c , where R c ss 10 km. 

2) Resonance layer (r — R Tes ). This is the 
region where the (energy-dependent) reso- 
nance condition in equation (4) is satisfied. 
When the adiabaticity condition 7 > 1 is ful- 
filled, this is also the region of emission of the 
sterile neutrinos. In our simulation we have 
studied the range R les /R c =0.8-1.8, namely 
we have assumed that the resonance layer is 
close to the core-atmosphere interface. 

3) Neutrino energy-sphere (r = Re)- The 
main reactions through which the v T ex- 
change energy with the stellar gas are NC 
scattering off electrons and neutrino pair 
processes. The neutrino energy sphere is 
defined as the region where the v T undergo 
the last inelastic interaction with the back- 
ground particles. For r > Re , neutrinos 
mainly scatter off nucleons. While this de- 
termines the neutrino transport opacity, the 
amount of energy exchanged is negligible, so 
that outside the energy-sphere the neutrinos 
start being thermally disconnected from the 
medium. 

4) Neutrino transport-sphere (neutrino sphere) 
(r — Rt)- This is the layer with optical 
depth « 1 which is determined by the neu- 
trino elastic cross section off nucleons. The 
cross section depends on the square of the 
neutrino energy, and thus Rt (as well as 
Re) is an energy-dependent quantity. Since 
Rt ^> Rics the initial asymmetry generated 
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at i? ros vanishes at Rt because of multiple 
scatterings that completely redistribute the 
momentum direction of the unconverted v r ■ 
For r > Rt the neutrinos no longer undergo 
diffusion processes and freely escape from 
the star. 

The density and temperature profile for the 
PNS core and atmosphere that have been used 
in our simulations are given by simple analyti- 
cal formulae. These expressions have been chosen 
according to the following criteria: (i) they ap- 
proximately reproduce the profiles resulting from 
detailed numerical studies of the PNS evolution 
(Burrows & Lattimer 1986); (ii) they satisfy the 
physical requirements that the v T decouple ener- 
getically with a spectral temperature T(Re) ~ 
8 MeV and around densities p(Re) ~ 10 12 g cm~ 3 ; 
(iii) they result in conservative estimates of the 
overall asymmetry in the sterile neutrino emission. 
The density profile is given by 



p(r) 




) (*r 



if r < R c ; 

if r > R c . 

(7) 

For suitable values of the exponents n c for the 
core and n a for the atmosphere, this expression 
reproduces reasonably well the results of a de- 
tailed numerical computation (Burrows & Lat- 
timer 1986), namely a slowly varying density for 
the core (n c ~ 1,2) and a steeper profile ~ r~ ria 
(with n a ~ 3-5) for the atmosphere. In Fig. 2 
we depict a set of profiles for different values of n c 
and n a . In our simulation we have used n c = 2 , 
which gives an appreciable density variation inside 
the core. With respect to a milder variation or to 
a constant profile, this is a conservative choice, 
since it favors the outwards diffusion of v r and re- 
duces the probability of random crossings of the 
adiabatic resonance region. For the atmosphere 
we used n a — 4 , which enhances with respect to 
steeper profiles, the probability that some of the 
neutrinos will diffuse back from the atmosphere 
and cross the resonance region in the 'wrong' di- 
rection. 

For the temperature profile we assumed a sim- 
ple model with an isothermal core 




T = 



T 



if r < R c ; 
if r > R c , 



(8) 



Fig. 1. — Variation of the adiabaticity parameter 
7 inside the resonance layer as a function of the 
angular distance from the magnetic dipole axis 
for different values of the decentred parameter 5 . 
Solid line: 6 = 0.1; dashed: 6 — 0.2; space-dotted: 
S = 0.3; dotted: S = 0.4; dash-dotted: <5 = 0.5 . 
The resonance is located at R rcs = 1.5 R c and the 
maximum value of the magnetic field strength at 
the resonance is B = 4 x 10 15 G. 
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Fig. 2. — Density profiles in units of 10 14 g cm' 3 
as a function of r/R c , for different values of the 
exponents n c and n a [see Eq. (7)]. The dark solid 
line depicts the profile that was used in the sim- 
ulations, and corresponds to n c = 2 and n a = 4 , 
the dashed line corresponds to n c = 1 and n a = 3 , 
and the light solid line corresponds to n c — 2 and 
n a = 5 . The position of the core-atmosphere in- 
terface (R c ) and the approximative positions of 
the v r energy-sphere (Re) and transport-sphere 
(Rt) are also shown. 
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with T = 30 MeV(Burrows & Lattimer 1986) and 
tit = 1. Since the ^ s are emitted with spectral 
temperatures close to the core temperature, larger 
values of T a would enhance the momentum asym- 
metry. During the early stage of the PNS forma- 
tion, the rapid deleptonization process produces 
changes in the temperature on the time scale of 
hundreds of milliseconds. In the first few seconds 
after core bounce, neutrinos are emitted predomi- 
nantly from the outer regions where the density is 
lower. Matter loses entropy, undergoes compres- 
sion, and the temperature increases to values much 
larger than the central temperature. The temper- 
ature inversion drives the heat flux towards the in- 
terior; however, the peak temperature reaches the 
central regions about 10 s after core bounce, when 
most of the neutrinos already escaped from the 
core (Burrows & Lattimer 1986). Since the neu- 
trino opacity grows with the square of the energy, 
a temperature inversion tends to keep the neutri- 
nos trapped in the inner regions, enhancing the 
probability for their conversion. We have verified 
that profiles with temperature inversion generally 
result in larger momentum asymmetries, and thus 
the isothermal core represents a conservative ap- 
proximation. As regards the temperature profile 
for the atmosphere, we assumed that it decreases 
with the radius as 1/r . With respect to steeper 
temperature gradients this again favors the prob- 
ability of 'wrong' back-crossings of the resonance. 
The D T elastic scattering off nucleons reads 

-^- = ^(C 1 + C 2 cose), (9) 
a cos v 8n 

where 9 is the scattering angle between the incom- 
ing and outgoing neutrino directions, and the coef- 
ficients C\ and C 2 depend on the neutrino-nucleon 
couplings averaged over the p and n number den- 
sities. For the PNS nuclear matter the canonical 
values C V A — —C\ — 1.26 given by isospin invari- 
ance are modified due to strange-quarks contribu- 
tion to the nucleon spin. We use the values C v = 



-1; 

c p A 



Cy 



1 - 4sin^6 w ^ 0.07, C r j 



-1.15. 



1.37 (Raffelt & Scckcl 1995; Keil, Janka, & 
Raffelt 1995). Using also Y p = 1 - Y n w 0.3 for 
the relative abundances of the nucleons, we obtain 

Ci=E ^v[(C£) 2 + 3(C£) 2 l«5.2,(10) 



Co 



N—n,p 

= E 

N—n,p 



Y, 



N 



-0.8. (11) 



Equation (9) determines the v T transport opacity. 
Inelastic reactions where the energy exchange is 
of order unity, such as neutrino-electron scatter- 
ing and neutrino-pair processes, give only a mi- 
nor contribution to the total opacity. However, 
they are responsible for keeping the neutrinos in 
thermal equilibrium with the background, and de- 
termine the relative positions of Rt and Re- We 
have assumed a relative rate between inelastic and 
elastic processes of the order of 10%, so once every 
ten scatterings we generate again the neutrino en- 
ergy according to the neutrino position within the 
star and to the corresponding local temperature. 
Since the v T emerge from the resonant layer with 
an overall momentum asymmetry equal and oppo- 
site to that of the v s , their thermalization and es- 
pecially the redistribution of their momentum due 
to multiple scattering is crucial for ensuring that 
the v s momentum asymmetry is left unbalanced. 

5. RESULTS 

Neutrinos are generated randomly inside the 
core with a Fermi-Dirac distribution correspond- 
ing to a spectral temperature T = 30MeV and 
zero chemical potential, and are left diffusing out- 
wards. Of course, the loss of v T due to conversion 
into sterile neutrinos eventually builds up a non- 
vanishing chemical potential. However, since our 
simulation of the neutrino diffusion is static (that 
is there is no time evolution of the relevant pa- 
rameters characterizing the PNS model) this ef- 
fect has been neglected. To study the effect of 
the decentred magnetic dipole, we have run sim- 
ulations with S ranging between 0.1 and 0.6. We 
have also studied the effect of moving the reso- 
nance layer from inside the core (R ICS /R C = 0.8, 
p ICS w 4.5 x p\4, vh Vb rs 50keV) to three different 
locations in the atmosphere R ICS /R C = 1.3,1.5 
and 1.8 (the last value corresponds to p ICS rj 
0.25 x pi4, m Vs w 12keV). The reference value 
for the neutrino magnetic moment has been fixed 
at p v = 10~ 12 pb (of course rescaling p v — > kp v 
and B — > B/k would leave the results unchanged) 
and we have varied the magnetic field in the star 
interior between 3 x 10 15 G and 8 x 10 15 G. The 
values of B correspond to maximum values of the 
adiabaticity parameter at the resonance surface 
in the range 7 w 3-12. In contrast, owing to 
the anisotropies of the magnetic field, in the res- 
onance regions far from the dipole centre, 7 <C 1. 
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Once the v T reach the resonance, the component 
of the magnetic field transverse to the direction of 
propagation is computed, and they are converted 
into v a with the appropriate transition probability 
Pv T ^v s - The v s escape from the star without in- 
teracting further. The large differences in the spa- 
tial values of 7 can result in huge anisotropics in 
the conversion efficiency between different points 
of the resonance layer, so that the v s are prefer- 
entially emitted from one hemisphere. In com- 
puting the v a momentum asymmetry we have ne- 
glected the effect of gravitational redshift. On 
the opposite hemisphere, most of the v T are not 
converted; however, further interactions with the 
medium wash out their asymmetry almost com- 
pletely. Of course, from regions close to the core- 
atmosphere interface, some v T can scatter inwards 
and resonantly convert when crossing the reso- 
nance region, so that a certain number of v a will 
be emitted in the 'wrong' direction, thus reducing 
the overall asymmetry. When the resonance lies 
inside the core, this effect is somewhat enhanced 
since a sizeable fraction of v T are produced in the 
region R ICS < r < R c . 

The fractional asymmetry in the total momen- 
tum of the neutrino emission needed to accelerate 
the remnant at a velocity V is given by 

^ * 03 ( v ) ( Etot ] 

Ptot ' U0 3 kms-V V3 x 10 53 crg S y \1.5M Q ) 

(12) 

where E tot is the total energy radiated in neutri- 
nos and M is the PNS mass. Assuming equipar- 
tition of the total luminosity between the active 
flavours, we see that to explain velocities of sev- 
eral hundreds of km s _1 the total asymmetry in 
the momentum carried by the v a and by the sur- 
viving T>t (p(v s + D T ) w 0.17ptot) should amount 
to w 10%-20%. The results of different simula- 
tions run with a statistics of w 4 x 10 5 v r are 
presented in Table I. It is apparent that momen- 
tum asymmetries of the correct size can be gen- 
erated. From the table it is also possible to infer 
the impact that the different parameters have on 
the size of the asymmetry. As expected, increas- 
ing the value of the decentred parameter decreases 
the fractional number of v s produced. The over- 
all asymmetry first increases; however, when the 
conversion efficiency decreases too much, also the 
asymmetry gets reduced. If the resonance is in- 
side the core, the asymmetries are generally small. 



This is mainly due to the large number of 'wrong' 
crossing of the resonance region, and to some ex- 
tent also to the moderate value of 7 max (larger 
values would require raising the magnetic field up 
to B rts 10 16 G). We stress that a crucial condi- 
tion that has to be satisfied to get large asymme- 
tries is that the adiabaticity parameter 7 should 
reach values sizeably larger than unity. In fact if 
7max ~ 1 the mechanism is rather inefficient, since 
only the D T with momentum almost orthogonal to 
the magnetic field direction will convert. Besides 
yielding a small conversion efficiency, this also im- 
plies that the momenta of the emerging v s will be 
mainly oriented tangcntially to the resonant layer 
and will not contribute to build up an asymmetry. 
Only for 7 max > 3-4, a sizeable momentum com- 
ponent in the radial direction can be generated. 

If the resonance is just outside the region where 
most of the neutrinos are produced (R ICS > R c ), 
asymmetries of the correct size are easily obtained. 
Because of the lower value of the local density, this 
allows for larger values of 7 max and also requires 
somewhat smaller magnetic fields. However, if 
[i v < 10 _12 ^s internal magnetic field strength 
B > 10 15 G are needed in any case. Since the 
relevant combination that determines the value of 
7 is the product /i„B, for larger transition mag- 
netic moments the required magnetic fields could 
, be smaller by a factor of a few. 

In spite of this, it is clear from our results that 
no obvious correlation between the kick velocity 
and the magnetic field strength can be expected. 
Similarly, any other correlation pattern would be 
quite difficult to recognize. In fact, for a given 
value of the neutrino transition magnetic moment 
/i„ the efficiency of the conversion is determined 
by a complicated interplay between the position of 
the resonance layer i? res , the value of the decentred 
parameter 5, and the strength of the magnetic field 
B. In some cases a large magnetic field together 
with a large value of S implies that 7> 1 only in 
a relatively small region, thus lowering too much 
the conversion efficiency and implying that the v s 
emission is mainly oriented tangcntially to the res- 
onance layer surface. Larger conversion rates are 
produced by smaller values of S, and sometimes 
this can also result in larger asymmetries. The 
exact value of the core temperature can also affect 
the results. For larger temperatures, about the 
same numbers of v s will be produced, but with a 



ff 



Table 1 
Results 



Rc 


B 


8 = ^e- 

Rc 


N(v T +v s ) 








0.8 


8 x 10 15 G 

(7max ~ 3.3) 


0.1 
0.2 
0.3 
0.5 


0.58 
0.43 
0.27 
0.08 


0.03 
0.06 
0.09 
0.14 


-0.001 
-0.002 
-0.004 
-0.004 


0.02 
0.05 
0.04 
0.03 


1.3 


6 x 10 lb G 

(7max ~ 5.5) 


0.2 
0.3 
0.4 
0.6 


0.57 
0.45 
0.33 
0.17 


0.11 
0.21 
0.29 
0.44 


-0.002 
-0.010 
-0.006 
-0.004 


0.09 
0.15 
0.17 
0.17 


1.5 


6 x 10 1S G 
( 7max ~ 11.3) 


0.3 
0.4 
0.5 
0.6 


0.51 
0.48 
0.44 
0.29 


0.13 
0.21 
0.27 
0.37 


+0.003 
-0.005 
-0.002 
-0.002 


0.10 
0.16 
0.16 
0.21 


1.5 


4 x 10 lb G 

(7max ~ 5.0) 


0.2 
0.3 
0.4 
0.5 


0.54 
0.44 
0.37 
0.26 


0.11 
0.19 
0.26 
0.36 


-0.005 
+0.001 
-0.006 
-0.003 


0.09 
0.12 
0.17 
0.15 


1.8 


3 x 10 1S G 

(7max ^ 7.0) 


0.2 
0.3 
0.4 
0.6 


0.60 
0.58 
0.43 
0.28 


0.09 
0.12 
0.22 
0.39 


-0.011 
+0.006 
+0.003 
+0.006 


0.07 
0.09 
0.15 
0.10 



Note. — The results of our simulations. Different positions of the resonance layer have been studied: 
inside the core (R TeB /R c = 0.8) and outside the core (R TeB /R c = 1.3,1.5,1.8). The values of the 
magnetic field are given in the second column, together with the corresponding maximum value of the 
adiabaticity parameter 7. In all the simulations we have assumed \i v = 10 — 12 [ib ■ The displacement of 
the magnetic dipole from the star centre <5 = Rg / R c is listed in the third column, while the resulting 
overall efficiency for P T conversion into u a is given in the fourth column. The momentum asymmetry 
for the v s , for the unconverted P T , and the overall v T + v s momentum asymmetry are given in the last 
three columns respectively. 
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higher spectral temperature. This effect tends to 
increase the asymmetry. 

Finally, we should mention that since the v a 
emission is likely to occur on a time scale much 
larger than the pulsar rotational period at birth, 
the resulting momentum kick would be averaged 
out proportionally to the cosine of the angle be- 
tween the spin axis and the magnetic axis. How- 
ever, even if we assume that the present orienta- 
tion of the magnetic axis is representative of its 
orientation at birth, we believe that it is not pos- 
sible to predict any unambiguous correlation be- 
tween the pulsars' velocities and their magnetic 
axis orientation, because of the large number of 
different parameters that concur to determine the 
overall effect. 

In conclusion, we have shown that v T <-» v s 
RSFP inside the PNS core is able to account for 
pulsar natal kicks of the required size. The cru- 
cial assumptions to achieve this result are (i) the 
presence of sizeable anisotropics in the PNS mag- 
netic field configuration; (ii) internal magnetic 
field strengths of the order of few x 10 15 G; (iii) a 
sterile neutrino with a mass of a few tens of keV 
and with a transition magnetic moment with an 
active flavour of the order of I0~ 12 fis- The first 
two assumptions have already been discussed in 
this paper, so let us comment briefly on the last 
point. If all the active neutrinos have a mass sat- 
isfying the cosmological limit < 100 eV, it would 
be somewhat complicated to construct a parti- 
cle physics model that realizes the third condi- 
tion, since a large magnetic moment generally im- 
plies rather large radiative contributions to the 
neutrino masses. This problem was extensively 
addressed in the past, and some clever solutions 
were proposed (Voloshin 1988; Barbieri & Moha- 
patra 1989). We believe that a consistent par- 
ticle physics model yielding a large v T -v s transi- 
tion magnetic moment and a light v T can be con- 
structed along the same lines. Let us also men- 
tion that quite recently the cosmological limits on 
neutrino masses have been revisited under the hy- 
pothesis that the Universe underwent a non-trivial 
thermal evolution right before the nucleosynthesis 
era (Giudice, Kolb & Riotto 2000). The result of 
this study opens up the possibility that also m VT 
could be of the order of a few tens of keV. Indeed 
this would render the whole picture much more 
natural from the particle physics point of view. 
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